Skin effect with arbitrary specularity in Maxwellian plasma 
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The problem of skin effect with arbitrary specularity in maxwellian plasma with specular-diffuse 
boundary conditions is solved. New analytical method is developed that makes it possible to obtain 
a solution up to an arbitrary degree of accuracy. The method is based on the idea of symmetric 
continuation of not only the electric field, but also electron distribution function. The solution is 
obtained in a form of von Neumann series. 
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I. INTRODUCTION 

The skin effect problem is one of the most important 
problems in plasma kinetic theory (see, for example, in 
Refs.[l[-@). The skin effect in plasma is a response of 
electron gas to external transverse electromagnetic field. 
The problem also has great practical importance. 

The solution of the skin effect problem with specular 
reflection boundary conditions is well-known [H , . The 
analytic solution of the problem with diffuse reflection 
boundary conditions has been obtained in the middle of 
the previous century (see, for example, Ref.Q). 

The skin effect problem with general specular - diffuse 
reflection boundary conditions [fj is not solved till now. 
It's well known that specularity coefficient q is a very 
important factor in the kinetic skin effect theory 0] - 
[§]. The limiting cases q = (diffuse surface scattering 
of electrons) and q — 1 (specular surface scattering of 
electrons) are only very special cases. Actually the spec- 
ularity coefficient q equals neither to zero, nor unit, and 
takes some intermediate values on. 

So, for example, in the work Q it is shown, that the 
specularity coefficient q is equal to 0.4 in Na wire. In 
this connection the skin effect problem with specular- 
diffuse boundary conditions has exclusively fundamental 
significance. Its value is great for theories, and for prac- 
tical applications. So it's obvious that the solution of the 
skin effect problem with general specular-diffuse reflec- 
tion boundary conditions is a very important task. 

The method of solution of this problem for degenerate 
plasma in metal has been developed in Ref. [8]. This 
method is based on the use of von Neumann series. Au- 
thors in Ref. [8] have demonstrated high efficiency of 
the method developed for computation of the skin effect 
characteristics. The goal of this work is generalization of 
the method developed in Ref. [§] in the case of gaseous 
plasma. 

By method of decomposition of the solution by eigen- 
functions of the corresponding characteristic equation ex- 
act solutions of the skin effect problem in metal for diffuse 



and specular boundary conditions are received in Refs. 

1, m 

In the last years interest to skin effect pro blems con- 
tinues to grow (see, for example, Refs. |ll|-[lrl|). In 
particular, in Ref. [lfj in limiting anomalous skin effect 
conditions, the oblique electromagnetic wave reflection 
from the sharp plasma boundary in an assumption of 
mixed (specular and diffuse) electron reflection from the 
boundary is considered. 



II. PROBLEM STATEMENT 

Let's gaseous (nondegenerate) plasma occupy a half- 
space x > . The distribution function of electrons 
/ = f(t, r, v) is normalized by electron numerical density 
(concentration of electrons): 



f(t,r,v) d 3 v = n(t, r), 



where p = tov is the electron momentum, m is the elec- 
tron mass, e is the electron charge, d 3 v = dv x dv y dv z . 

We consider electromagnetic wave which propagates 
in direction orthogonal to the plasma surface. Then the 
external field has only one y -component. The internal 
field inside plasma has only y -component 



Ey(t, X ) 



*E(x) 



too, where u> is the field frequency. 

To describe the electron distribution function we will 
use the Vlasov — Boltzmann kinetic equation. The colli- 
sion integral will be represented in the form of r -model 

df df _, . . df fo(v) — f(x,v x , t) 
at ox op y t 

where r is the time between two electrons collisions, r = 
v is the effective electron collision frequency, fo(v) 
is the equilibrium maxwellian distribution function, 



2^ 



k R T 



3/2 



exp 



mv s 
2fc B T, 
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For weak fields this equation may be linearized: 

/ = f (C) (l + C y h(x, C x )e-^<) . (1.1) 

Here 

/ Fi\ 3/2 

k{C )=n(L) cxp(-C 2 ), 

where (3 = m/ (2k B T) , C = y/fiv , C is the dimension- 
less electron velocity, fcs is the Boltzmann costant, T 
is the plasma temperature. 

For function h(x,C x ) we have the following kinetic 
equation: 

^ dh(xu») +ZQh[xull) = e{x ^ (L2) 



where 



z = l- i 



In the equation (1.2) x\ is the dimensionless coordi- 
nate, 

H = v xV /f j3, xi = j = v v r f3x, tl=Vt, 

ti is the dimensionless time, / is mean electron free path 
and e(.Ti) is the dimensionless electric field: 

e(xi) = W^T E(Xl) - 

We will neglect displacement current. Then the equa- 
tion for electric field may be written in the form: 



e = — —3y( Xl )> 



(1.3) 



where j y {xi) is the electric current, 

Jv(xi) = e J v v [f (v) + y/pv y e- iut h(x,v x )}d 3 v. (1.4) 

Let's extend the electric field and the electric distri- 
bution function on the "negative" half-space x < in 
symmetric manner. For the functions E(x) and h(x,^i) 
we will then write: 

e(xi) = e(-xi), h{xi,fi) = h{-xi,-^). (1.5) 

We may rewrite the equation (1.4) with use of dimen- 
sional parameter a : 



d 2 e(xi) . a 
dx 2 

Here 



= -i—= [ cxp(— [i' 2 )h(xi, [i') dy! . (1. 
l V n J 



6) 



47re 2 n£j „/l\ 2 



a = 



c 2 [3v 3 m \5 



S 2 = 



, l = V T T, V T = — 

mv 



1 



vt is the electron thermal velocity, 5 is the penetration 
length of external electric field for normal skin effect, a 
is the plasma conductivity. 

By extension procedure (1.5) on the half-space x < 
we may include the surface conditions in the equation for 
skin effect problem. 

Specular - diffuse boundary conditions on the bound- 
aries of positive and negative half-spaces may be written 
in the form: 



h(+0,fj,) = qh(+0,-fj,), < fi < 1, 



/i(-0,/i) = qh(-0,-fj,), -1 < n < 0, 

where q is the specularity coefficient, < q < 1 . 
In accordance with (1.5) we obtain: 

h(+0, n) = qh(-0, n), < fj, < 1, (1.7) 



h{-0, /i) = qh(+0, /i), -1< n < 0. (1.8) 

The required function h(xi,fi) and the electric field 
must decay away from the surface: 



h(+oo,n) = 0, e(+oo) = 0. 



(1.9) 



We assume that the gradient of the electric field is 
finite and known at the plasma boundary: 



e'(0) = 4 |e;|<+oo. 



(1.10) 



Here, the gradient of the electric field on the plasma 
boundary e' s is given. 



III. CHARACTERISTIC SYSTEM 

The variable xi will be denoted again by x . 

Let's include boundary conditions (1.7) and (1.8) in 
the kinetic equation (1.2), and boundary condition (1.10) 
include in the electric field equation (1.6). 

As a result we will obtain system of equation for skin 
effect in half-space of the plasma: 

dh , . . 
11 dx + z ° h ^ x ^> = 

e(x)-(l-q)\n\h(T0,»)5(x), ±/x > 0, (2.1) 



2-KLoa ' 



d 2 e{x) 
dx 2 



CO 

= — i~7= / cxp(— [i' 2 )h(x, fi') d/j,'+ 
V n J 
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2e s 5(x). (2.2) 
The impedance is determined by formula [![: 

A-Kiu E y (0) 



Z = 



c 2 E'joy 



With the use of dimensionless field e(x) this relation 
may be rewritten in the form: 

_ Amujl e(0) 

Z = j ~- 

c e s 

From equation (2.1) and boundary conditions (1.9) we 
obtain the following expression for x > 0, /i < : 

oo 

h + (x, a) = cxp f — ") / exp f — ^ e(t) dA. 

A« v At 1 J v A* ' 

a: 

In the case x < 0, /i > we obtain: 

X 

h- (x, fx) — — exp ( — ) / exp ( — ) e(t) dt. 

At \ A 4 ' J v A* ' 

— oo 

Then we may rewrite the equation (2.1) in the form: 

9h , , , , , 
ii— + 2ott(^, M) - e(x) = 

-(l-q)\lJ,\h ± (0,fi)S(x), ±/i>0. (2.3) 

The solution of the equations (2.2) and (2.3) we will 
seek in the form of Fourier integrals: 



oo 

' ) = h j elkXE ^ dk > 



(2.4) 



2tt 



exp(ikx)E{k) 
z + ik/j 



dk. 



(2.6) 



It's may be proved, that the expression for h-{x,ji) 
coincides with the expression for h + (x,/i). Therefore 
we have 



MO, AO 



oo 

i r E{k)dk 



2n J zq + ik/i 



We substitute the expressions (2.4), (2.5) and (2.6) 
into the equations (2.2) and (2.3). This procedure leads 
to characteristic system of equations: 

$(fc, n)(zo + iky) = 



*(*)-(!-< / 



\y\ f £(fci)d/ci 



Zq + iklfj, ' 



(2.7) 



k^E(k) = 2e s - ia / exp(-f/)^(k,fi)dti. (2.. 



The function e(x) is an even function. Then E(— k) 
E(k) , and equation (2.7) may be rewritten as 



$(&,/*) (so +ikn) 



E(k)-(l-qy- 



fj,\zo f E{k\)dk\ 



rr J Zq ■ k'/r 




(2.9) 



Let's substitute the expression (2.9) into the equation 
(2.8). Then we obtain: 



<5(x) = — / e™*dk, 



L{k)E{k) = -2e s + 



h(x,fi) = — I e tkx <I>(k lf j,)dk. 
2tt 



(2.5) 



Then for the function h+(x,n) the following expres- 
sion may be derived: 

2n(j, 



(1-q)— t / E(h) J{k,ki)dkx. (2.10) 



Here L(k) is the dispersion function, 



2 .a f exp(— jjL 2 ) dji 



L(k) = k 2 -i 



IT J Zq + ikfi 



oo 

dkjdt^{ikt 



z t 



E(k) 



k 2 - 



2izoa f cxp(— /j, 2 ) d/j, 
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and J(k,ki) is the integral 



J(k,h) = -?=/" 

J ( z o 



exp(— u 2 ) du 



+ k 2 u)(z 2 + k 2 uY 



Characteristic system consists of two equations (2.9) 
and (2.10). 

The integral J{k\, k 2 ) we will express through the in- 
tegral exponential function. We will spread out integrand 
partial fractions: 



1 



A B 

+ 



(zq + k\u){zl + k\u) Zq + k\u z 2 + k 2 u' 



where 



A = — 



k\ 



B= k * 



7 2/p_L2\' " ~iii.i_i.iy 
Thus, we receive, that 



du 



B 







or 



where 



J(k u k 2 ) = AJ (h) + BJ (k 2 ), 



2 

& cx p(^2 ) y ex p(- ^2 )— . j = !. 2 - 



IV. SOLUTION OF THE PROBLEM WITH THE 
USE OF VON NEUMANN SERIES 

Let's expand the solution of equations (2.9), (2.10) by 
the following series: 

E(k) = E {k) + (l-q)E 1 (k) + (l-q) 2 E 2 (k) + --- , (3.1) 



${k,[i) = & {k,ii) + (1 - <?)$i(fc,M)+ 



(l-g) 2 $ 2 (fc, M ) 



(3.2) 



Functions Ej(k) and ^j(k,fi) (j = 1,2,3,- ••) may 
be obtained from the characteristic system. For zero ap- 
proximation we have: 



L(ky 



L(k)(z + ikjj) ' 



(3.3) 



For first approximation we obtain: 



Ei(k) 



L(k)-K r i 



E a (ki)J(k,ki)dki 



(3.4) 



Mm 



E (ki)dk 



zo + ikfi (z + ikfi)ir J z 2 + k\\i 2 
o 



\. (3.5) 



For n -approximation the following expression may be 
derived: 

E n (k) = 



n r 

J E n _ 1 (k 1 )J(k,k 1 )dk 1 , n=l,2,-", (3.6) 
E n (k) 



L(k) 



z + ifc/u 



(20 + i 



00 

Ml f E^xjk^dkx 

ikfi)-K J Zq + k\p? 



n = l,2,- 



(3.7) 



We may rewrite the expressions (3.1) in the form 



' az 2 /" J(fc, fci) 



L(k)L{k x ) 



E 2 {k) = -2e(^)\ 

\ 7TZ / 

J{k,k 1 )j{k ll k 2 ) 



dk\, 







L(fc) J L(fci)L(fc 2 ) 



dkidk 2 



In general case when n = 1, 2, 3, • • • , we have: 
£?„(*) = -2e s (— ^ x 

J{k,k\)J{k\,k 2 ) ■ ■ ■ J(k n -i,k n ) 



00 00 







L(fc)L(fci) • • • L(fc„) 



rffci • • • 



(3.8) 

Therefore the series (3.1) constructed may be ex- 
pressed in the explicit form: 



E(k) = 



2e c 



n=l 



00 00 







L(k) 

J{k,ki)J(ki,k 2 ) ■ ■ ■ J(k n -i,k n ) 
Hh) ■ ■ ■ L(k n ) 



dki ■ ■ ■ dk„ 
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V. ELECTRIC FIELD, DISTRIBUTION 
FUNCTION AND SURFACE IMPEDANCE 

In accordance with (2.4) and (2.5) we will construct 
expressions for electric field and distribution function. 
Using the expressions (3.1) and (3.2) we obtain: 

oo 

i °° r 

e(x) = -^2(1- q) n / E n (k)coskxdk, (4.1) 



c e s J 
o 

We decompose Z in the following series: 
Z = Z + (1- q)Z 1 + (1 - q) 2 Z 2 + ■ ■ 

Here 

oo 

4 '"' ' E n (k)dk, n = 0,1,2,' 



(4.4) 



(4.5) 



i r 

h(x,v) = -Y / (l-qT e ikx $ n (k,[i)dk. (4.2) 



We rewrite the expression (4.1) with the use of (3. 
in the following form: 



Now we write down expressions for zero, first and sec- 
ond approximations for impedance (see (4.4) and (4.5)) 



Z = 



AiluI f dk 8iwl f dk 



L(k) 



o 



(4.6) 



, . 2e c f coskxdk 





z, = — 



iu)l az 2 . f dk f J(k,ki)dk 



r- ,T I J L(k) J I. ( ; ) 





oo oo 



J(k 1 k 1 )J{k u k 2 ) ■ ■ ■ J(fc„_i,fc„) 
L{k!) ■ ■ ■ L(k n ) 



dk\ ■ ■ ■ dk n 



o o 

Function h{x,jj) may be written in the form: 



h(x,fx) 



1 

27 



E(k)- 



Zo\fJ,\ f E(ki) dk\ 



,T J Zq + lyi 1 





e lkx dk 



(4.3) 



z + ikfx 

Expression (4.3) may be written also in the next form: 

oo 

1 f zq coskx + kfismkx 



h(x,n) 



2q + k 2 [i 2 



zo|a*| f E(ki) dki 



E(k)- 



IT 2 J Zq + k\[l 2 



dk. 



If we know function h(x, fi) we may write down the 
electron distribution function / according to the equa- 
tion (1.1). 

Let us consider now the calculation of impedance: 



z = We(0) = 2^Z i 



8iul razQ\ 2 



c 2 \ 



TV I 



dk f J(k,ki)dki f J(ki,k2)dk 



L(k) J L(fci) J L(k 2 ) 
o o 

Expression for general term of series (4.4) has the form: 

8iusl / azg \ " 
it i J 



c 2 \ 



oo oo 







J( fc,fcl)J(fci,fc 2 ) ■■ ■ J(fc ra -l,fcn) 

L(k)L{h) ■ ■ ■ L{k n ) 
x dk dk\ ■ ■ ■ dk n ■ 

VI. ANALYSIS AND DISCUSSION 

In the previous sections we have considered the 
method, leading to exact solution of the skin effect prob- 
lem with arbitrary specularity coefficient. In the case 
q = 1 the method leads to the classical solution (4.6) 
of the problem with specular surface conditions (see, for 
example [H, [lfj], [13] )• ^ n this classical solution is 
represented in the form: 



Z t 



8iujI 



■of 



d,T 



\(izor) ' 



(5.1a) 
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where 



X(iz r) 



a t 



cxp(— ji 2 ) d[i 

/i — IZ{)T 



(5.16) 



The comparison of the expressions (5.1) and (4.6) gives 
Z re f = Zq . Indeed, after the change of variables in the 
integral we have: 



dr 



\{iz T) 



dr 



1 A(«2 /i 



dk 





Now let us consider second approximation: 

Z = Z + (l-q)Z 1 + (l-q) 2 Z 2 . 

When q = 1 this solution is exact. Maximum devi- 
ation from exact solution corresponds to the case when 
q = . The exact solution of the problem in the case 
q = is also well known [Io| : 



Zdif - 



A'Klbjl 



r OO 






\{iz Q T) 


L 





dr 

7? 



(5.2) 



It's convenient to rewrite the expression (5.2) with the 
use of our notations: 



Z dif 



In 



k- 2 L(k) 



dk 



The curves 3 on the plots correspond the ratios 
of impedance for diffuse scattering surface condition 
to impedance for specular scattering surface condition 
M(Z dif )/$t(Z ref ) (%(Z dif )/$(Z ref ) ). 

The derived method has maximum error in the case of 
extremely anomalous skin effect, when parameter a 3> 
1 . In this case ratios defined above are equal to 1.125. So 
it is obvious from the plots, that in zero approximation 
the method error is equal to 12.5% . 

For the first approximation in this case we have 
Z d if I Z = 1-03 • So the first approximation error is equal 
to 3% . For the second approximation in this case we 
have Z^f I Z = 1-01. And the second approximation 
error is equal to 1% . 

The analysis of plots shows, that the considered 
impedance ratios in first approximation coincide with the 
exact solution when a < 10 _1 . For the second approxi- 
mation the coincidence is observed when a < 1 . 



VII. CONCLUSION 

The effective method of the solution of boundary prob- 
lems of the kinetic theory is developed. This method is 
based on symmetric continuation of the electric field and 
distribution function of electrons. 

The offered method gives an error less, than 1% in the 
second approximation already. The method is accurate 
and allows to construct exact solution in the form of von 
Neumann series. 



Now consider the ratio of real (and imaginary) parts of 
the solutions constructed in zero, first and second approx- 
imations to the solution in zero approximation ^R(Zq) 
( ^s(Zq) ) for the case q = . The last solution coincide 
with the solution of the problem with specular scattering 
boundary conditions. 

We will build two plots (curves 1 and 2): 



M(Z ) 



= 1 



g(gi) 
ft(Z ) 



and 



or 



Y 2 = 



jft(Z + Z x + gg) 



= 1 



jft(Zx) 9K(Z 2 ) 

n(z ) n(z ) 



Y 2 =Y 1 



$t(Z 2 ) 

n{z Q ) ' 



and also analogous plots for the ratios of the imaginary 
parts. 

Here Zq is the solution for the case of specular suface 
conditions. Values Z\ and Z 2 correspond to corrections 
for the first and the second approximations. 



Re(z)/Re{z ) 




FIG. 1: Dependence of 3?(Z)/3J(Zo) on parameter a for the 
case q — . The curve 3 corresponds to diffuse scattering 
boundary conditions, the curves 1,2 correspond to first and 
second approximations. 
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